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I Introduction 

This report is the final report for the Cooperative Agreement NCC-2- 
1000. The tasks outlined in the various proposals are listed below 
with a brief comment as to the research performed. The publications 
are attached. 


II. NCC 2 1000 Tasks and Accomplishments 


1. Development of an aerosol chemistry model (Tasks 1-5, Dr’s 
Tabazedeh and Lin) 

This was the main focus of this Cooperative Agreement and the 
results are contained in two published papers attached) 

2 Utilization of satellite measurements of trace gases along with 

analysis of temperatures and dynamic conditions to understand ice 
cloud formation, dehydration and sedimentation in the winter polar 
regions. (Task 6 Dr Stone) 

This was a specific task for Dr. Stone and resulted in a paper has been 
submitted for publication (attached) 

3. Comparison of the HALOE and SAGE II time dependencies of the 
Pinatubo aerosol decay (Task 7, Dr. Bergstrom) 

This was a specific task for Dr. Bergstrom and resulted in a paper that 
is being prepared for publication. 
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Abstract. Recent laboratory data on ice freezing from 
aqueous H 2 S0 4 solutions are used to update our homogeneous 
ice freezing nucleation algorithm. The effect of particle size on 
the ice freezing process is demonstrated by calculating and 
comparing ice freezing curves for different size particles 
ranging from 0.2 to 4 micron in radius. Using the ice 
nucleation model we show that liquid water saturation is 
required above -44°C to activate submicron H 2 S0 4 particles 
into cloud droplets. A thermodynamic model is used to show 
that the available laboratory data on ice freezing from ternary 
H 2 S0 4 /HN0 3 /H 2 0 solutions are insufficient to adequately 
address the effect of HN0 3 uptake on the homogeneous freezing 
of aqueous H 2 S0 4 panicles into ice. 

Introduction 

Ice clouds play important roles in Earth’s climate 
[Ramanathan el ai , 1983] and chemistry Solomon [1990]. 
One mechanism that has been suggested to initiate ice cloud 
nucleation in the atmosphere is homogeneous freezing of 
aqueous H^SC^ solution droplets and Dodd, 1989; 

DeMoti el al., 1994; Jensen el al., 1994; Heymsfield and 
Miloshevich , 1995; Demon el al., 1997]. However, the ice 
freezing process through this mechanism is complicated by 
the fact that a significant fraction of gas phase HN0 3 will 
partition into sulfate aerosols near the ice freezing point of the 
solution [Tabazadeh el al., 1994; Carslaw el al., 1994]. 

The ice freezing curve in aqueous H 2 S0 4 droplets was 
recently measured by Koop ei al. [1998]. The ice freezing 
temperatures determined by Bertram ei al. [1996] for aqueous 
H 2 S0 4 droplets, which we previously used to parameterize our 
ice nucleation model [Tabazadeh el al., 1997 a, b], appear to be 
systematically too high [Koop et ai, 1998]. In addition, 
differential scanning calorimetry has recently been used to 
determine the ice freezing curve in aqueous HN0 3 and ternary 
H 2 S0 4 /HN0 3 /H 2 0 droplets [Chang ei al . , 1999], In this work 
we use recent laboratory measurements on ice freezing from 
binary and ternary solution droplets [Koop et al., 1998; Chang 
ei al., 1999], along with new surface tension data on the 
aqueous H 2 S0 4 system [Myhre ei al., 1998], to update our 
homogeneous ice freezing nucleation algorithm. 
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Model Update 

Figure 1 shows critical ice saturations expected to nucleate 
ice in aqueous H 2 S0 4 , HN0 3 and ternary H 2 S0 4 /HN0 3 /H 2 0 
solutions as a function of the ambient water vapor pressure 
based on recent laboratory measurements. The data of Koop ei 
al. [1998] and Bertram ei al. [1996] show very different 
behaviors for the freezing of ice from aqueous H 2 S0 4 solution 
droplets. The main objective of this paper is to update our ice 
nucleation model to agree with more recent measurements. 

An interesting feature of Figure 1 is that the ice freezing 
curves for both aqueous H 2 S0 4 and HN0 3 solution droplets are 
nearly identical. According to equilibrium model calculations 
[Tabazadeh ei al . , 1994; Carslaw ei ai, 1994], the ionic 
strengths of H 2 S0 4 /H 2 0 and HN0 3 /H 2 0 at the same relative 


HjSO^ (Koop et al.) - - - H^SO^ (Bertram et al.) 

HN0 3 (Chang et al.) - - Nucleation Model 



Figure 1. The variation of critical ice saturation as a function 
of water partial pressure over the solution. The solid and dotted 
curves are calculated directly from the tabulated relations given 
in Koop et ai [1998] and Chang et ai [1999]. The numbers on 
the dotted lines give the H 2 S0 4 content of a ternary 
HN0 3 /H 2 S0 4 /H 2 0 aerosol in weight percent The solid dotted 
curve is based on the data of Bertram et ai [1996] convened to 
an appropriate form using the method described in Tabazadeh 
et al. [1997a]. Also shown are the results of our ice nucleation 
model for a 4 micron aqueous H 2 S0 4 particle. 
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humidity arc similar in magnitude. Since the equilibrium 
freezing point depression is often directly proportional to the 
ionic strength, the fact that the two binary electrolytes have 
similar ice freezing kinetics is at least partially rationalized 

Also shown in Figure 1 are the ice freezing curves for 
H : S0 4 /H NO ; /H : 0 ternary solutions, which contain I and 5 r r 
H : S0 4 by weight [C/uwjr a al .. 19991. The slight i< 2W) 
increase in the critical icc saturation in a ternary solution, as 
compared to the binary systems, is rationalized by the increase 
in entropy of the solution as a result of mixing. Below we use a 
thermodynamic model to show that the current laboratory data 
on icc freezing from ternary solutions are insufficient to 
adequately address the effect of HNO* uptake {by sulfate 
aerosols) on the ice formation process in the atmosphere. 

Denoting L, >cm^) and as the volume and weight percent 
of a H : S0 4 solution droplet, the rate of ice nuclealion (particle 
sec' 1 ) at a given temperature is: 


molecules in solution based on laboratory see freezing data 
We follow the same procedure here to derive a new relation for 
the activation free energy of water molecules in solution hased 
on the ice freezing data of Koop w ah 1 1998], 

In deriving an expression for the activation free energy, wc 
generated two fittings of the activation energy function for 
particle sizes of 1.55 and 6.3 micron in radius that bracket the 
laboratory size range [ Koop a ah. 1998], The final activation 
energy relation reported in Table 1 is obtained by averaging 
the two expressions above for different particle size 
assumptions. In addition, we replaced the surface tension 
relation in the model by a new expression based on recent 
laboratory measurements by Myhrc ct al. |I998]. All the 
functions required for the nuclealion calculations are described 
and summarized in Table 1. The nuclealion model calculates a 
unique solution for the ice freezing temperature of a given size 
droplet as a function of only the ambient w ater vapor pressure. 


J = G T. 




F k {T, w\ s )- A F„ (/ (7\uy ) 

-- — 


(1) 


where C (particle sec 1 ) is the preexponential factor, A F 
(ergs) is the Gibbs free energy for the formation of the ice 
germ, &F (H , tergs) is the diffusion activation energy of water 
molecules across the icc/sulfale solution phase boundary, and 
k is the Bo!i2mann constant. The preexponenlial factor and 
the Gibbs free energy are estimated as follows: 

C(T. . V d ) = 2. 1 x 10” V d > w/lV >; v ,T)7- (2) 


4 > 

a g,' 


(3) 


where <7 ulf , ut . is the interface energy, between the ice/sul face 
solution and (cm) is the critical germ radius. For our 
purposes, an fee freezing event occurs when J = I particle sec' 1 . 
We previously described [Tabazadeh e\ al., 1997 a, b] how 
the equations (1) through (3) can be used to extract a numerical 
relation for the variation of activation free energy of water 


Atmospheric Implications 

Figure 2 shows the effect of particle size on ice nuclealion 
from an aqueous H : SCE solution droplet. The size-dependence 
calculated by our model is in good agreement with the size- 
dependence measured in the laboratory for ice nuclealion in 
pure supercooled water droplets (sec Figure 2 in Pmppachcr 
[1995]). As expected, due to smaller volumes ( Vj term in 
equation 1), higher ice saturations are needed to form ice in 
submicron atmospheric panicles as compared to micron-sized 
aerosols studied in the laboratory [Koop et al.. 1998]. Recent 
measurements [Chen ei ah. 2000] support our calculations 
.that submicron sulfate panicles freeze at lower temperatures as 
compared to micron-sized particles. 

Laboratory freezing data on submicron water droplets 
indicate that such droplets can often remain supercooled to 
temperatures near and sometimes below 1 -40 C [Hagen er ah. 
198); Pmppacher . 1995]. Thus the fact that our ice nuclealion 
model calculates a liquid water activation temperature below 
-40 "C for submicron aqueous H : S0 4 particle> (Figure 2) is not 
surprising since the effect of solute is to further depress the 
freezing point of ice in solution beyond that of a pure 
supercooled water droplet. Overall, a nuclealion model offers 


Table 1. Ice Nuclealion Functions 


Lateni heat of melting in erg mol 1 (7* is in K): L ni (T) - 10 7 [6005.2356+ 1 8.27 19t T -273. 1 5) - 0 . 06354 : 7- 273 15)* j 
lee density [T is in K > . p u ,AT) = 0.916 - 8.75 * \ 0~*(T - 273.15) - 1 667 x I0“"(7 - 273.15)* 

Sulfate soluiion/air interface energy based on the data of My lire el al. [1998] in dyn cm' 1 (where w is the H ; S0 4 weight percent of the solution): 
aJ*J J , /j|r = 85 755071 14 + 9.54196631 8.v 1 0“ 2 iv - 1.103647657x1 0" 1 h * + 7.485866933xl0~ 1 n ? - 1 .912224 1 54.vl0~ J w 4 + 1 736"S9787.\10^iv 5 
o~" air = 82.01 1 97792 + 5.3 1 2072092a 10 _1 w - 1. 050692 1 23*10" V + 5.415260617 a I0" ? h 1 - 1 . 1455"3827.vl0" J u 4 + 8.969:57061.vl0' 7 w 5 
= 77.40682664 - 6.963 123274.vlO -? w - 9.682499074.vl0~ ? u" + 8.87979880 a 1 0 1 tv* - 2.3846695 16x10"' n- 4 + 2.09535S048.vl0" 7 if 5 
Fur 1 80 < T < 220 K. a T ul , Mr = a;,f lmr + (5.5 - 0.025n<<;L, - ) 

For 220 < T < 260 K. o] u)lai , = <7™.,,, + (6.5-0.025r>(a^)„ (r - 0 %)„ ir ) 

Diffusion activation energy x I0 n in ergs derived from the Koop ei at. [1998] ice nuclealion data (7' is in K): 

Fur T < 220 K. iF,„ = -17459.5 1 6 1 83 + 458.4582755 IT - 4.849283 1 3 I IT 2 + 0.0260036588787' 5 - * : 99 1577798 x1 O' 5 C 
+8.90490946 1 R.vl O'* C - 2 .49322574 ] 9.vl 0' 1 1 T* 

For T > 220 K = 104525. 93058-1103. 764465IT+ l.0703327027" : + 0 01 73862543227' •' - I 5506854268 ilO - '"/" 1 

-3 ,266 1 9 1 24974:10"* T* + 6 4679.54459.vl 0~ 111 T* 

Equilibrium sulfuric acid weight percent composition and water saturation vapor pressure: See Taba:a%u h ei al. [1997a] 

Sulfuric acid soluiion density: See Myhrc et al. [1998] 
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Figure 2. The variation of the critical ice saturation (and 
supercooling) as a function of water partial pressure and 
particle size in radius calculated using the ice nucleation 
model. The critical supercooling is defined as the difference 
between the critical ice freezing temperature and the ice frost 
point of the solution. Liquid water activation temperatures 
marked on the plot show transition temperatures above which 
particles activate into liquid water droplets prior to freezing 
into ice. The. variable X in the above equation is In (PH : Oi. 
where PH ; 0 is the ambient water vapor pressure in mb. 


several ad\antages compared to parameterization* for treating 
the ice formation process in the atmosphere. Simple 
parameterization^ tend to provide information that is 
exclusive to one size particle, whereas in the atmosphere ice 
nucleation in principle occurs over an entire submicron 
aerosol particle size distribution. 

Liquid water activation temperatures marked in Figure 2 i for 
three different particle sizes* show transition temperatures 
above which sulfate particles will first activate into liquid 
water droplets prior to freezing into ice. To accurately predict 
this transition point in the atmosphere is important since 
cirrus ice clouds that fiisi go through a liquid water activation 
phase are often composed of a large population of small ice 
crystals [./cu.vc// cf ni. Radiative properties of cirrus 

clouds are strongly affected b> the number of ice particles 
nucleated \Jcn<cn and huai, 1994], The results of our 
calculations show that submicron ILSt.V particles in the 
atmosphere can remain supercooled to temperatures near -44 
C Therefore it is likely that above -44 C cirrus ice clouds 
olicn form by first going through a liquid water activation 
phase provided that they are devoid of heterogeneous solid 
impurities \Marnn. 1998; Dttnoii el ai, 199"). 


The composition of an aqueous H : S0 4 solution droplet in the 
atmosphere is altered by the uptake of HNO v In Figure 3 the 
change in the ternary solution composition as a Junction ot 
temperature is shown for three selected water vapor pressure 
profiles, which bracket typical conditions in the lower 
stratosphere and upper troposphere. The shaded area shows the 
range of H : S0 4 compositions of ternary systems for w hich the 
ice freezing properties of solutions have been determined in 
the laboratory [Chang el ai , 1999]. 

Ternary solution droplets containing between I to 5 wt % 
H : S0 4 freeze at ice saturations near that of the binary systems 
[Chang ei ai, 1999, see Figure 1]. In the background 
stratosphere the H^S0 4 weight percent in solution near the ice 
nucleation point is located within the gray area for which the 
ice freezing properties of ternary solutions have been studied 
in the laboratory. Thus the effect of HN0 3 uptake in depressing 
the ice freezing temperature is negligible (Figure I j if the 
stratosphere is in a background state (Figure 3 k However, for a 
volcanically perturbed state, where H : S0 4 mixing ratios are 
higher, the H : S0 4 weight percent of a ternary solution near the 
ice freezing point is close to 15 (Figure 3). Because this 


Background H 2 S0 4 wt % - - Perturbed H ? S0 4 wt % 

Background HNC^ wt % - • Perturbed HNG 3 wt % 



Temperature (K) 


Figure 3. The variation of ternary H : S0 4 /HN0 ; /H : 0 aerosol 
composition with temperature for different water vapor 
pressure regimes based on the thermody namic relations given 
in Lnn ei ai [1995]. The 0 1 and 0 01 mb curves are calculated 
for a constant H : S0 4 atmospheric mixing ratio of |()0 ppt and 
variable HNO ; atmospheric mixing ratios of 0.1 (background) 
and 2 ppb (periurbed by convective pollution) at 200 mb. The 
0,0003 mb curves are calculated for a constant UNO, 
atmospheric mixing ratio of 10 ppb and variable H : S0 4 
mixing ratios of 0.5 (background) and 20 ppb {perturbed by 
volcanic eruptions) at 60 mb. The calculations were performed 
down to aboul -3.2 K below the equilibrium ice condensation 
point of t he solution. The shaded area shows i he range of 
H : SO, compositions in ternary aerosol solutions for which the 
ice freezing properties of ternary sweats have been measured 
in the laboratory | Chan a et ai . I9vq. The arrows m the plot 
refer to H_S0 4 composition in ternary solutions for which 
laboratory ice freezing data are lacking. Laboratory freezing 
data tor ternary solutions are currently unavailable for volcanic 
and background conditions in the stratosphere and upper 
t ro pos p he re , respect i vc I y 
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ternary composition regime has not been studied in the 
laboratory, it is difficult to compare and contrast the effect of 
HN0 3 uptake on the ice freezing process in the stratosphere 
between the two cases. Satellite data indicate that it may be 
more difficult to freeze aerosols in a volcanically perturbed 
environment as compared to a background state [Santee et al . , 
1998], Thus additional laboratory data are needed to better 
understand how volcanic eruptions may affect the ice 
formation process in the lower stratosphere. 

Similarly in the upper troposphere the laboratory data on ice 
freezing form ternary systems are insufficient to adequately 
describe and compare the differences between the background 
and perturbed states. For example, in a polluted upper 
troposphere HNO ; mixing ratios reach values as high as 2 
ppbv [Laaksonen et al . , 1997; Schnieder el al . , 1998]. The 
high levels of HNO ; in the upper troposphere are linked to dry 
convection that injects polluted boundary layer air directly 
into the upper troposphere. The H 2 S0 4 weight percent of 
ternary solutions in a polluted upper troposphere are located 
within the gray shaded area (see Figure 3). On the other hand, 
ice freezing curves for a background state (between 8 to 12 % 
H 2 S0 4 and I to 5 % HNO 3 by weight), containing about 100 
pptv of HNO :i , have not yet been explored in the laboratory 
(Figure 3). Thus it is difficult to address the question of how 
convection (leading to high levels of HN0 3 ) may affect cirrus 
formation in the upper troposphere when no information is 
available on ice freezing properties for compositions 
representative of the background state (-100 pptv of HN0 3 ). 

Summary 

We updated our homogeneous ice freezing nucieation code 
using recent laboratory measurements. Our results indicate that 
ice freezing curves for submicron aerosol particles are 
significantly different than those determined in the laboratory 
for micron-sized particles. Further we show that in the 
atmosphere, for temperatures above -44 C, submicron aqueous 
H : SOj particles will activate into liquid water droplets prior to 
freezing into ice. Finally the available laboratory data on icc 
freezing from H : S0 4 /H N0 3 /H 2 0 solutions are insufficient to 
address how HNCu uptake by aqueous H 2 S0 4 particles may 
affect the ice formation process in the atmosphere. 
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Abstract. Simple expressions are fitted to the results obtained from ion interaction 
thermodynamic models for calculating HNO3 and H2O vapor pressures over 
the NH3/H2SO4/HNO3/H2O system at cold temperatures. The vapor pressure 
expressions are incorporated into a mass conserving equilibrium solver for computing 
aerosol compositions in the lower stratosphere and upper troposphere. The 
compositions calculated from the aerosol physical chemistry model (APCM) are 
compared against previous parameterizations. The APCM compositions are in 
better agreement with the compositions obtained from ion interaction models than 
from other previous formulations of the NH3/H2SO4/HNO3/H2O system. The only 
advantage of the APCM over the ion interaction approach is that the numerical 
scheme used in the model is fast and efficient for incorporation into large-scale 
models. The APCM is used to calculate HNO3 solubility in ammoniated aerosols 
as a function of HNO3. H2SO4 and NH3 mass loadings in the lower stratosphere 
and upper troposphere. While the uptake of HNO3 by ammoniated aerosols is 
strongly dependent upon the solution neutrality (or pH), we find that in both 
the lower stratosphere and upper troposphere a significant fraction of HNO3 will 
exist in aerosol solutions near and below the ice frost point irrespective of solution 
neutrality. 


1. Introduction 

Thermodynamic electrolyte models are often used for 
calculating properties of inorganic aerosols in the lower 
troposphere [Stetson and Seinfeld. 1981; Pilinis and Se- 
infeld \ 1987; Wexler and Seinfeld , 1991; Kim and Se- 
infeld , 1995; Jacobson et a/., 1996; Jacobson, 1999b] 
and stratosphere [ Tabazadeh et ai. 1994: Carslaw et ai, 
1994, 1995b; Weisenstein et ai , 1997]. In the lower tro- 
posphere. aerosol models are often used in air quality 
studies to assess the effects of aerosols on health, gas- 
phase partitioning, and visibility. In the stratosphere, 
aerosol models have been used to simulate the forma- 
tion and growth of polar stratospheric clouds, which are 
linked to stratospheric ozone depletion [ Solomon , 1999]. 
However, thermodynamic treatments in large-scale at- 
mospheric models are in general not suited for calcu- 
lating aerosol compositions in the upper troposphere. 
Since upper tropospheric aerosols participate in the nu- 


cleation and growth of cirrus clouds and may also be 
involved in the scavenging of trace gas species, it is 
important to understand their chemical and physical 
properties [Karcher and Solomon , 1999]. 

Recently, equilibrium aerosol formulations have been 
incorporated into three-dimensional models to simulate 
the radiative impacts of aerosols on climate [ Adams et 
ai , 1999; Jacobson , 2000]. In this work, an aerosol phys- 
ical chemistry model (APCM) with an efficient solv- 
ing scheme suitable for incorporation into large-scale 
atmospheric models is developed. The APCM is com- 
pared against various parameterizations, including re- 
sults from the aerosol inorganics model (AIM2) of Clegg 
et ai [1998a]. 

2. Background on Thermodynamic 
Aerosol Models 

2.1. Ion Interaction Approach 

The ion interaction approach is originally based on 
the work of Pitzer [1991]. The solution behavior in 
Pitzer’s method is determined by a series of single- 
ion (and water) activity equations that are based on 
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Table 1. A n Coefficients for the Vapor Pressure of HN0 3 a 


r 

ao 

ai 

a 2 

a 3 


a 5 

a 6 

a 7 

0.0 

12.6923 

- 0.0325979 

- 0.071719 

0.812331 

2.34482 

0.000438486 

- 0.00428908 

- 0.000709075 

0. 1 

12.-4691 

- 0 . 00 1 9955 

- 0.077484 

0.727085 

2.42858 

0.000259723 

- 0.00442994 

- 0.000706633 

0.2 

12,2326 

0.0260933 

- 0.085355 

0 . 64-1593 

2.52519 

9 . 32848e-05 

- 0.00453571 

- 0.000694942 

0.3 

11.9840 

0.05 13596 

- 0.094561 

0.567113 

2.63058 

- 5 . 53402e-05 

- 0.00461276 

- 0.000678339 

0.-4 

1 1.7205 

0.0725417 

- 0.104727 

0.300403 

2.74278 

- 0.000172282 

- 0.00466394 

- 0.000659102 

0.5 

1 1.4289 

0.0550553 

- 0.116820 

0.450799 

2.86783 

- 0.000241534 

- 0.00468708 

- 0.000632214 

0.6 

11.0837 

0.0953196 

- 0.134271 

0.415462 

3.02666 

- 0.000267273 

- 0.00467490 

- 0.000579027 

0.7 

10.6697 

0.106923 

- 0.160517 

0.378195 

3.24022 

- 0.000284624 

- 0.00461980 

- 0.000480710 

0.3 

10.2087 

0 . 11-295 

- 0.194708 

0.324242 

3.50279 

- 0.000327867 

- 0.00451884 

- 0.000342092 

0.9 

9.73753 

0.130202 

- 0.232997 

0.251054 

3.79040 

- 0.000406493 

- 0.00437409 

- O . O0O184569 

1.0 

9.28445 

0.144365 

- 0.271608 

0.165018 

4.07962 

- 0.000510598 

- 0.00418932 

- 2 . 90206e-05 

1.1 

8.87335 

0.156600 

- 0.306440 

0.080649 

4.34517 

- 0.000611458 

- 0.00396949 

0.000101871 

1.2 

8.55964 

0 . 15"392 

- 0.327751 

0.039034 

4.52665 

- 0.000618013 

- 0.00372273 

0.000153819 

1.3 

8.33474 

0.134427 

- 0.331971 

0.092212 

4.60197 

- 0.000411043 

- 0.00345255 

0.000106943 

1.-4 

7.95881 

0.110127 

- 0.354342 

0.142888 

4.78919 

- 0.000208202 

- 0.00313268 

0.000158658 

1.5 

7.51749 

0.0956133 

- 0.390173 

0.143658 

5.05815 

- 0.000120908 

- 0.00275735 

0.000282545 

1.6 

7.14929 

0.0560305 

- 0.421553 

0.118501 

5.29796 

- 0.000104343 

- 0.00234767 

0.000379673 

1.7 

6.88502 

0.0761817 

- 0.441786 

0.092775 

5.46818 

- 0.000107787 

- 0.00192436 

0.000413991 

1.8 

6.68262 

0.0630763 

♦ 0.453078 

0.080818 

5.58436 

- 0.000101314 

- 0.00149773 

0.000397780 

1.9 

6.43124 

0.0455338 

- 0.465952 

0.086457 

5.71507 

- 7 . 48198e-05 

- 0.00106043 

0.000386668 

2.0 

h 

5.72588 

0.0242624 

- 0.524014 

0.096528 

6.13955 

- 5 . 31470e-05 

- 0.00053951 

0.000609824 


1 1.5869 

0.222271 

- 0.109389 

0.279027 

2.83237 

- 0.00101574 

- 0.00488557 

- 0.000704950 


a A n - ao + aiwj - a 2 w 2 + a 3v /^7 + a 4y /w£ + a 5 wf + a 6 wiw 2 + a 7 w|. 
b The oc represents the ternary system of NH 3 /HN0 3 /H 2 0 (i.e., H 2 S0 4 = 0). 


thermodynamic properties of mixed solutions. Physi- 
cal parameters for the fundamental activity relations in 
an ion interaction model are fitted to laboratory mea- 
surements of solute-water mixtures of interest. Ther- 
modynamic electrolyte models based on the ion in- 
teraction approach have evolved significantly in recent 
years by Clegg and coworkers [e.g.. Clegg and Brim- 
blecombe, 1990, J995a. 1995b: Carslaw et ai , 1995a; 
Clegg et ai. 1998a, 1998b]. Since the models of Clegg 
and coworkers, collectively referred to as aerosol in- 
organics model (-AIM2), are fitted to laboratory mea- 
surements conducted over a broad temperature range 
(< 200 to 328°K), they are more accurate than those 
based on the common binary activity approach de- 
scribed in section 2.2. However, the complex nature 
of the ion interaction approach makes AIM2 computa- 
tionally impractical for three-dimensional applications. 
Here we use AIM2 to generate solution compositions 
and vapor pressures for the NH 3 /H 2 S0 4 /HX0 3 /H>0 
system for a wide range of humidities and temperatures. 

2.2. Binary Activity Approach 

The second approach first emerged from modeling 
inorganic aerosols in air quality studies (e.g.. Stetson 
and Seinfeld. 1981; Bassett and Seinfeld , 1983; Saxena 
et ai , 19SGL This method separates out water (via 


either the Gibbs-Duhem equation or the water equa- 
tion [Stokes and Robinson, 1966]) and solute activity 
coefficients. Usually, a mixing rule of either Bromley 
[1973] or Kusik and Meissner [1978] is used to estimate 
the mixed solute activity coefficients. The equilibrium 
models based on this approach require only a knowl- 
edge of water and solute activity coefficients at a binarv 
level. Thermodynamic models based on the binary ac- 
tivity approach are computationally more efficient but 
less accurate than the ion interaction approach because 
the physics of the activity coefficients in the former 
approach are mainly based on the behavior of binarv 
solutions instead of mixed solution properties used in 
the latter formulations. To contrast two different ther- 
modynamic treatments, we will update an equilibrium 
model based on the binary activity approach, EQU- 
SOLV II [Jacobson et ai , 1996; Jacobson , 1999b], and 
compare its predictions against results obtained from 
AIM2. EQUISOLV II applies a well-converged numeri- 
cal solver scheme to simultaneously solve a large number 
of equilibrium equations. 

2.3. Vapor Pressure Approach 

A combination of the two approaches outlined in sec- 
tions 2.1 and 2.2 is used in this work to develop a 
fast and accurate parameterization of the NH 3 /H 2 SO^/ 
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Table 2 . B n Coefficients for the Vapor Pressure of HNOs a 


r 

bo 

b . 

b 2 

b 3 

b < 

b 5 

b6 

b 7 

0.0 

- 7672.46 

67. "280 

32.3655 

- 259.196 

- 235.30 

- 0.256884 

0.97336 

0.197376 

0.1 

- 7609.77 

50.7068 

32.9809 

- 222.190 

- 255.43 

- 0.131368 

1.07699 

0.201805 

0.2 

- 7541.22 

35.7444 

34.4801 

- 189.780 

- 280.54 

- 0.026085 

1.16451 

0.202182 

0.3 

- 7467.24 

22. "371 

36.6718 

- 162.151 

- 309.62 

0.059118 

1.23812 

0.199536 

0.4 

- 7387.11 

11.9765 

39.4988 

- 140.155 

- 342.40 

0.122744 

1.29896 

0.194190 

0.5 

- 7297.89 

3.4748 

43.2297 

- 124.533 

- 380.56 

0.163341 

1.34688 

0.184702 

0.6 

- 7194.29 

- 3.0516 

48.6478 

- 113.773 

- 428.87 

0.184446 

1.38064 

0.166777 

0.7 

- 7074.27 

- 8.6193 

56.4190 

- 103.509 

- 491.39 

0.196255 

1.39907 

0.136740 

0.8 

- 6944.62 

- 14.0870 

66.1301 

- 90.235 

- 565.46 

0.207582 

1.40222 

0.096867 

0.9 

- 6815.18 

- 19.5349 

76.6808 

- 73.739 

- 644.27 

0.220387 

1.39134 

0.053245 

1.0 

- 6693.25 

- 24.6807 

87.0206 

- 56.230 

- 721.33 

0.231211 

1.36800 

0.011651 

1.1 

- 6584.97 

- 28.3719 

96.0022 

- 42.304 

- 789.58 

0.230911 

1.33398 

- 0.021582 

1.2 

- 6504.53 

- 27.8717 

100.995 

- 44.004 

- 832.71 

0.192747 

1.29170 

- 0.031795 

1.3 

- 6449.42 

- 19.3426 

100.908 

- 77.488 

- 843.95 

0.079582 

1.24188 

- 0.012861 

1.4 

- 6347.58 

- 8.4904 

106.038 

- 117.723 

- 886.93 

- 0.052980 

1.17584 

- 0.022354 

1.5 

- 6203.79 

- 0.1957 

118.040 

- 143.506 

- 971.43 

- 0.156090 

1.08786 

- 0.069346 

1.6 

- 6055.69 

4.3272 

133.047 

- 152.440 

- 1073.17 

- 0.221191 

0.98027 

- 0.132468 

1.7 

- 5927.18 

6.5545 

148.101 

- 146.787 

- 1173.72 

- 0.249961 

0.85834 

- 0.195765 

1,8 

- 5829.77 

6.2290 

161.504 

- 129.145 

- 1262.57 

- 0.245338 

0.72779 

- 0.250386 

1,9 

- 5768.38 

3.1411 

172.326 

- 101.381 

- 1334.05 

- 0.209351 

0.59425 

- 0.291726 

2.0 

- 5744.66 

- 2.5250 

180.115 

- 63.965 

- 1385.49 

- 0.141562 

0.46318 

- 0.317848 

oc b 

- 7298.14 

- 51.7642 

35.6176 

- 55.070 

- 360.093 

0.376735 

1.41665 

0.252861 


a B n = bo + bj wj - b 2 w 2 + b 3 ^/wf + b4 ywj 4- bswj + w 2 H- b?w^, 
b The cc represents the ternary system of NH3/HNO3/H2O (i.e., H2SO4 = 0). 


HNO3/H2O system for atmospheric applications. The 
vapor ^pressure approach requires generation of large 
sets of H 2 0 and HNO3 vapor pressure arrays over the 
NH3/H2SO4/ HNO3/H2O system for many assumed 
compositions and temperatures. The vapor pressure, 
composition, and temperature fields generated are fitted 
into simple mathematical expressions similar to those 
derived by Luo et al. [1995] for the H2SO4/HNO3/H2O 
system. The vapor pressure relations are then in- 
sorted into a numerical mass-conserving equilibrium 
solver similar to that used in EQITSOLY II [Jacobson 
tt a /.. 199G] for gas-aerosol phase partitioning calcula- 
tions. 

3. Model Development 

3 . 1 . Aerosol Physical Chemistry Model 
( APCM) 

Luo et al. [1995] have shown that equilibrium par- 
tial pressures of HNO3 and H 2 0 over the H2SO4/HNO3/ 
HA) system roughly follow a Clausius-Clapeyron rela- 
tion of the form InP = A + B/T (where P is pressure. 
T is temperature, and A and B are constants) for a 
fixed solution composition. For the quaternary system 
of NH3/H2SO4/HNO3/H2O we found a similar behav- 


ior for the variation of HNO3 and H 2 0 vapor pressures 
over the solution. 

Assuming that NH 3 and H2SO4 reside completely in 
the condensed phase (i.e., both NH 3(g ) and H^SO^g) 
are negligible), we follow the approach of Luo et al. 
[1995] and define wi and w 2 as the weight percents of 
ammoniated sulfate and nitric acid, respectively, in the 
solution as follows: 

W] = weight % of (NH4) r H2-rS04(aq), 

w 2 = weight % of HN0 3{aq) , (1) 

NH 3 = NH ^q) (2) 

r H 2 SO, HSO- aq) +SO^ aq) ’ 

where r in w t is the fixed mole ratio of ammonia to sul- 
furic acid, and (aq) is the aqueous-phase species. The r 
ratio can be considered as the degree of the ammoniated 
solution neutrality and can take on any values in APCM 
(including fractions) between 0.0 and 2.0. For example, 
if r equals 0.0, 1.0, and 2.0, w L represents the weight 
percents of H2SO4, (NH^HSO-j, and (NH^SOj, re- 
spectively, in the solution. Using Cleyy et alJ s [1998a] 
model, we first generated a series of arrays of HN0 3 and 
H 2 0 vapor pressures for a wide range of weight percent 
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Table 3. C w Coefficients for the Vapor Pressure of H 2 O a 


r 

Co 

Cl 

C2 

c 3 

c 4 

c 5 

c 6 

Cr 

0.0 

22.7490 

0.0424817 

0.0533280 

- 0.0567432 

- 0.276555 

- 0.000621533 

- 0.000311769 

- 0.000283120 

0.1 

22.7669 

0.0310031 

0.0493720 

- 0.0270045 

- 0.265333 

- 0.000496930 

- 0.000261912 

- 0.000257919 

0.2 

22.7838 

0.0223810 

0.0457014 

- 0.0059721 

- 0.254223 

- 0.000403349 

- 0.000224673 

- 0.000234421 

0.3 

22.8005 

0.0165207 

0.0423183 

0.0065020 

- 0.243407 

- 0.000339520 

- 0.000198990 

- 0.000212594 

0.4 

22.8179 

0.013221 5 

0.0393297 

0.0110504 

- 0.233647 

- 0.000303073 

- 0.000184277 

- 0.000192973 

0.5 

22.8367 

0.0119853 

0.0369883 

0.0096303 

- 0.226435 

- 0.000288721 

- 0.000180339 

- 0.000177085 

0.6 

22.8556 

0.0119008 

0.0354244 

0.0058633 

- 0.222617 

- 0.000287316 

- 0.000186497 

- 0.000165709 

0.7 

22.8706 

0.0122336 

0.0342757 

0.0024719 

- 0.220187 

- 0.000291401 

- 0.000201103 

- 0.000156776 

0.8 

22.8791 

0.0129460 

0.0329580 

- 0.0009337 

- 0.215893 

- 0.000300122 

- 0.000221710 

- 0.000146949 

0.9 

22.8818 

0.0143893 

0.0311831 

- 0.0063234 

- 0.208231 

- 0.000316451 

- 0.000246179 

- 0.000134615 

1.0 

22.8812 

0.0168823 

0.0289977 

- 0.0153842 

- 0.197574 

- 0.000343217 

- 0.000273061 

- 0.000120122 

1.1 

22.8803 

0.0204917 

0.0267376 

- 0.0284821 

- 0.185921 

- 0.000381016 

- 0.000302055 

- 0.000105483 

1.2 

22.8846 

0.0243030 

0.0255578 

- 0.0417119 

- 0.179964 

- 0.000421023 

- 0.000333949 

- 9 . 75631 e -05 

1.3 

22.8946 

0.0259818 

0.0267000 

- 0.0449953 

- 0.186824 

- 0.000440711 

- 0.000370519 

- 0.000103697 

1.4 

22.8950 

0.02651 82 

0.0275800 

- 0.0428170 

- 0.191549 

- 0.000449334 

- 0.000409584 

- 0.000108446 

1.5 

22.8866 

0. 0275456 

0.0272887 

- 0.0426774 

- 0.189011 

- 0.000462731 

- 0.000447757 

- 0.000106656 

1.6 

22.8752 

0.0293967 

0.0263037 

- 0.0461528 

- 0.182140 

- 0.000484748 

- 0.000483688 

- 0.000101242 

1.7 

22.8640 

0.0319034 

0.0251308 

- 0.0524783 

- 0.173920 

- 0.000514332 

- 0.000517346 

- 9 . 51958 e -05 

1.8 

22.8542 

0.0347733 

0.0240901 

- 0.0603014 

- 0.166228 

- 0.000549028 

- 0.000549196 

- 9.0361 le -05 

1.9 

22.8459 

0.0376949 

0.0233346 

- 0.0682171 

- 0.159987 

- 0.000586063 

- 0.000579624 

- 8 . 75512 e -05 

2.0 

22.8381 

0.0403666 

0.0228641 

- 0.0749269 

- 0.155242 

- 0.000622692 

- 0.000608700 

- 8 . 66960 e -05 

DC b 

22.7279 

0.0163541 

0.0303799 

- 0.0208336 

- 0.151376 

- 0.000354967 

- 0.000653371 

- 0.000197299 


a C w = Co + Cl Wi + C 2 W 2 + C 3v /w7 + C4>/W2 + CsW^ + CeW^ + C 7 W.f. 
b The DC represents the ternary system of NH3/HNO3/H2O (i.e., H2SO4 = 0). 


combinations (w h w 2 ; where 0.0% < wj + w 2 < 85.0%) 
at cold temperatures, ranging from 180 to 270°K. The 
vapor pressures were then fitted into the following func- 
tions: 


0.0% < w j w* < 85.0%, 180°K < T < 270°K, (3) 


in v / , , B„(r l w ll w 2 ) /r 

1uPhno 3 = A n (r. wi,w-)) + ^ , (4j 


IhPh.o = C w (r. w \ . w 2 ) + 


T 

Pw(r.wi,w 2 ) 

T 


(0) 


where Phno 3 and Ph^o are in inbar. Coefficients A„. 
B n . C w . and D w were then parameterized into a mathe- 
matical function of wj and w 2 similar to that described 
by Luo el al. [1995]. The fitting results are summarized 
in Tables 1-4 for different r ratios, ranging from 0.0 and 
2.0 (r values are separated by ari increment of 0.1). 

In Plate 1 the calculated vapor pressures using equa- 
tions (4) and (5) are compared with those obtained from 
AIM2 for three r ratios (shown in Plates la, lb, and 1c). 
In general, Ph 2 o calculated by equation (5) agrees well 
with the model of Clegg et al [1998a]. However, the 


agreement of Phno 3 using equation (4) is rather poor 
for most compositions. As shown in Plate 1, the cal- 
culated Phno 3 becomes increasingly inaccurate as the 
solution neutrality (i.e., r) increases. For some cases the 
difference between estimated Phno 3 values and those 
predicted by AIM2 reached as high as a factor of 5. To 
ensure differences are minimized, Diff terms account- 
ing for the differences (deviations) of the “calculated” 
values (estimated from equations (4) and (5)) and the 
“real” values (obtained from AIM2) were added to equa.- 
tions (4) and (5): 

10 \ t \ , B„(r,w h w 3 ) i 

InP HNO3 = A„(r, wt, w 2 ) + + 


Diff„(r, w [ . w >, T) 

in o , \ , Dw(r,w,,w*>) , 

lnPuaO = C w (r, w h w 2 ) + j + 

Diff w (r, w j . \v 2 , T). 


( 6 ) 

(7) 


For a fixed i\ W|, and w 2 the Diff correction terms ap- 
proximately follow simple polynomial functions in tem- 
perature: 

Diff,, (r, \vi , w 2 , T) = cr n0 + a nl T + cr n2 T*’ (8) 
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Table 4. D w Coefficients for the Vapor Pressure of H- 2 O a 


r 

do 

di 

d 2 

d 3 

d<i 

d 5 

d 6 

d 7 

0.0 

-5850.24 

21.9744 

1.48745 

-44.5210 

59.6400 

-0.384362 

-0.644671 

-0.208556 

0.1 

-5849.94 

19.8594 

1.73322 

-39.7129 

58.0581 

-0.336545 

-0.606959 

-0.210361 

0.2 

-5849.15 

17.5776 

1.92716 

-34.3704 

56.5456 

-0.289409 

-0.568784 

-0.211815 

0.3 

-5848.04 

15.1511 

2.06677 

-28.5548 

55.1548 

-0.243024 

-0.530295 

-0.212915 

0.4 

-5846.79 

12.6262 

2.12829 

-22.4377 

54.0535 

-0.197733 

-0.491539 

-0.213544 

0.5 

-5845.48 

10.1120 

2.06175 

-16.4657 

53.5369 

-0. 154503 

-0.452501 

-0.213397 

0.6 

-5843.72 

7.79449 

1.85460 

-11.3889 

53.6970 

-0.115030 

-0.413300 

-0.212396 

0.7 

-5840.63 

5.79437 

1.60810 

-7.6468 

53.9687 

-0.0803938 

-0.374279 

-0.211123 

0.8 

-5835.71 

4.05810 

1.46205 

-4.9098 

53.5707 

-0.0500282 

-0.335935 

-0.210378 

0.9 

-5829.33 

2.44674 

1.47151 

-2.5108 

52.2170 

-0.0225318 

-0.298645 

-0.210471 

1.0 

-5822.28 

0.84270 

1.60302 

0.0831 

50.1179 

0.0033213 

-0.262577 

-0.211194 

1.1 

-5815.44 

-0.79036 

1.74377 

2.9989 

47.9281 

0.0280087 

-0.227587 

-0.211883 

1.2 

-5810.31 

•2.20859 

1.55496 

5.1396 

47.5949 

0.0492844 

-0.193207 

-0 2 10535 

1.3 

-5807.19 

-2.71615 

0.62592 

3.4557 

51.4599 

0.0605252 

-0.158668 

-0.204729 

1.4 

-5801.66 

-2.43122 

-0.39261 

-1.5150 

55.7388 

0.0628630 

-0.124306 

-0.198357 

1 .5 

-5792.68 

- 1 .85557 

-1.06065 

-7.4818 

57.9126 

0.0612095 

-0.091134 

-0.193964 

1.6 

-5781.23 

-1.28487 

-1.31402 

-13.0971 

57.6578 

0.0586697 

-0.059876 

-0.191880 

1.7 

-5768.33 

-0.88563 

-1.18923 

-17.6257 

55.2142 

0.0571507 

-0.030919 

-0.191856 

1.8 

-5754.83 

-0.75269 

-0.74364 

-20.6765 

50.9331 

0.0578653 

-0.004421 

-0.193536 

1.9 

-5741.30 

-0.93290 

-0.03072 

-22.0791 

45.1436 

0.0615518 

0.019549 

-0.196612 

2.0 

-5728.09 

-1.44363 

0.91523 

-21.7913 

38.0640 

0.0686466 

0.040959 

-0.200889 

CC b 

-5736.70 

-0.27462 

2.02241 

0.5073 

20.6699 

0.0337051 

-0.045872 

-0.179320 


a D w = d 0 + d i w i + d 2 w 2 + d 3v Av7 + d 4 ywj + d 5 w’f + d 6 wiw 2 + d 7 wf. 
b The oc represents the ternary system of NH3/HXO3/H2O (i.e., H2SO4 = 0). 


- Diff w (r. wj . wo. T) = a w Q + <r wl T + <7 w 2 Tn (9) 

where (7 are second order polynomial coefficients fitted 
for a particular combination of wj, w>. and r. The 
effects of including Diff terms in improving the vapor 
pressure fits are shown in Plates Id- If where HNO3 and 
HoO vapor pressures calculated from equations (6) and 
(7) (with polynomial coefficients of equations (8) and 
(9)) are compared against AIM2 results for the same r 
ratios examined above. As shown in Plates Id- If, Diff 
terms force the calculated vapor pressures of both H 2 0 
and HN0 3 to agree with AIM2 results. For r values 
other than 0.0. 1.0. and 2.0, including Diff terms in 
vapor pressure relations produces a nearly exact agree- 
ment with AIM2 results. Thus, for APCM we tabulated 
polynomial coefficients of a (lookup tables available as 
electronic supporting material 1 ) for all possible combi- 
nations of w t . w 2 . and r. The weight percents of w t and 


'Supporting lookup tobies are available via Web browser or via Anonymous 
FTP from ftp: Vvwmos.jgu.org. direelor> r “apenJ ,, (L'scmanie v "anonymous \ 
Password ® "guest”): subdirectories in the (tp site are arranged by paper 
number Information on searching and submitting electronic supplements is 
found at http wvvu.agu org pubs esupp_about.html. 


w 2 in the lookup tables cover and span the composition 
spectrum from 0.0% to 85.0% with 1% increments in 
composition. 

As temperatures cool, solubilities of trace gases, such 
as nitric acid, increase significantly, thereby depleting 
gas phase concentrations. Thus it is essential to cal- 
culate the distribution of HNO3 between the gas and 
aerosol phases. To simulate the gas-aerosol partition- 
ing of HNO3, equations (6) and (7) are coupled with 
mass conservation of HNO3: 

total HNOj = Phnoj + HN0 3(aq ) (in mol/ra* 1 ). (10) 

Three unknowns, (NH.t) r H 2 _ r SO. Uaf|) (or w t ), HN03 (aq ) 
(or w 2 ), and Phno 3 - can be uniquely determined by 
equations (6), (7), and (10). For this work, the same 
numerical scheme as that utilized in EQl'ISOLY II [ Ja- 
cobson et al. } 1996] is applied to solve equations (6)-( 10) 
iteratively. 

At iteration steps where W| and w 2 are not tabu- 
lated in the lookup tables, the bilinear (area weighted) 
averaging [Jacobson. 1999a] is adopted for interpola- 
tion. For example, to estimate properties at [wj , w 2 ] = 
[25.3. 10.8], results obtained from four adjacent points, 
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Table 5. Composition Functions a 




A 

B 

c 

D 




(NH A ) 2 S0 4 


0.00 < a w 

< 0.30 





yi 


5. 1229560233e+01 

-9.6193408336e-01 

1.472057870ie+02 

-1.72291 i2367e-02 

y> 


4.5462384945e+01 

-9.48720289 I7e-01 

3.34121551660+02 

-2.1 t24Cl4825e+02 

0.30 < a w 

< 0.60 





yt 


9.734 59 16855e — 01 

-2.7444251275e+00 

-l. 89665286 14c+01 

1 .92 1 2056647e+0 1 

\’2 


1.0l4S647694e+00 

-2.39 1696 1554e+00 

-2.29879509860+01 

2. 1920018865e+0l 

0.60 < a w 

< 1.00 





yt 


-5.37563894 90e-t- 02 

1 .3556456453e-02 

-1.8951 24828 le+01 

5.5692 105238e+02 

Y‘2 


-4. 2921 43l877c+02 

1.2875856340e-02 

-2.1275933520e+01 

4.5075478762e+02 




NH a NOs 


0.00 < a w 

< 0.40 





yi 


1 .31 65236 l33e+02 

-1. 0103415600c 4-00 

-3.l726829968e+0l 

-1.40420757610+02 

ya 


8.1 149051028e+01 

-9.6213385243e-0l 

1.0536766694e+02 

-1.5360919227C+02 

0.40 < a w 

< 0.85 





yt 


8.5782004669e+0l 

-I.2795635979e+00 

2.99826614420+01 

-1.0776290885e+02 

>2 


3.2723371 194e+01 

-1,32992178510+00 

-8.12459730750+00 

-2.45202617990+01 

0.85 < a w 

< 1.00 





yt 


—1.496959951 6e — 01 

-2.0271382540c+0l 

-2. 1065828892c+02 

2.0679076661e+02 

y > 


1.66360621870+03 

-1.3^65428608e— 01 

1.8592287444e+02 

-1.84966766350+03 




nh a hso 4 


0.00 < a w 

< 0.25 





yi 


8.914S556445e+0l 

-9.2225449765e-01 

2. 1668 15771 2e+ 02 

-2.75523385310+02 

\'2 


5.98408745 14c+01 

-9.3729948048e-01 

2.72262245220+02 

-2. 1 72 1584622e+02 

0.25 < a w 

< 0.80 





yi 


2.8756459725e+00 

-2.55071 56 l09e+00 

-3.33 1 6860933e — 0 1 

- 1. 849421 40 12e— 03 

ya 


3.5550024880e+00 

-2.09103505000+00 

-1.04925398290+01 

7.89230135136+00 

0.80 < a w 

< 1.00 





yi 


- l.4227339985e+00 

2.9277231432e+01 

- 1 . 72 1 5058088e + 0 1 

l.8441662636e+01 

y-j 


-7.0483234 774e-0l 

2.0832745270e+01 

-2.3289593929e+01 

2.3901 137559e+0l 


a v = Aa®. + Ca w + D Read 5.1229560233e+01 as 5. 1 229560233 x 10 l . All parameterization are 
valid for 190° K < T < 260° K only. Molality m is calculated by m(a w ,T) = y i (a w ) + (T — 190) 
[(/■j(a u ) - y i (a a -)] /TO. where a w is water activity (relative humidity expressed in fraction). Compo- 
sition functions for HjSCU and HXO3 in the same mathematical form are given by Tabazadeh et al. 
[1997a, 19.97 bj. 


[25.0. 10.0]. [25.0. 11.0], [26.0. 10.0], [26.0, 11.0], are 
weighted and averaged, depending on how far the de- 
sired point is located away from the four fixed points. 
Similarly, a simple linear interpolation (of averaging 
predictions from two adjacent r ratios) is implemented 
if initial r falls inberween the 0.1 increments (see Ta- 
ble 1). Since the intervals are small (1% for \v t and w-> 
and 0.1 for r), the two interpolation methods yield rea- 
sonable results that are in close agreement with AIM2 
in predicting HXO :J and HjO vapor pressures over the 
XHj/HjSO i/HXOj/HjO system. 

For conditions where LRSCb mixing ratio is zero we 
define w t and \v-> as the weight percents of NH.jX 03 | afj j 
and HXOj ( j UI }. respectively, and repeat the above steps 
to obtain vapor pressures and solution compositions for 
the system. 


3.2. An Update of EQUISOLV II 

In addition to APCM parameterization an equilib- 
rium model based on the binary activity approach [Stet- 
son and Seinfeld , 1981; Pilinis and Seinfeld, 1987; Wexi- 
er and Seinfeld , 1991; Kim and Seinfeld , 1995; Jacob- 
son et a/., 1996; Jacobson , 1999b], EQl ISOLV II [Ja- 
cobson et al . , 1996; Jacobson , 1999b], is updated for 
application at colder temperatures. The updated EQ- 
UISOLV II will be compared against AIM2 and APCM. 
The important equilibrium equations to solve in EQUI- 
SOLV II for the NH 3 /H.>SO,/HNO;,/H,.0 system are 


HSO 


m«t} 




HNO-i 


3(g) 


u + 


H r«n, + xo 


* 


(ii) 


( 12 ) 
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Table 6. Mean Binary Activity Coefficients** 


<*i(T) = 00 + 0 iT + &r 2 

Coefficients 

Pa 

Px 

Pa 




H NOz{n — b) 




1.1550607531e+0l 

-8.6598065340e- 02 

1.6 100335 132e-04 

Ot[ 


-2.6045646534e+01 

1. 85701181 52e-01 

-3.403293 164 le-04 

a 2 


9.6192566768e+00 

-6.4314120728e-02 

1.1363959743e-04 

Q 3 


-1.4330104905e+00 

9.3470490393e-03 

— 1 .632 1 073928e— 05 

a 4 


9.518176781 8e— 02 

-6.1 14013U74e-04 

1. 06058569 19e-06 



-2.3568483 179e-03 

M936021604e-05 

-2.5755879088e-08 




{NH<) 7 SO A {n=Z) 




3.44283 16946e+00 

-4.5328321654e-02 

1 .0202 137529e-04 

a 1 


-l.0030630804e+00 

6.0339372495e-03 

- 1 .2 162329645e-05 

q 2 


4.2925175339e — 02 

-3.0262230896e-04 

5.9904 176768e-07 



-6.68l3265363e-04 

5.003 1688798e-06 

-9. 74 6025858 7e-09 




NH<NO z (n=4) 


a 0 


2.7037375056e+00 

- 1 . 7880 1 7 1 4 98e-02 

2.8497778393e-05 



-4. 942468584 4e+00 

2. l621849334e-02 

-2.2789057453e-05 

a 2 


9.0584578176e— 01 

-4 .294 7 1 6 1 383e-03 

4.7971905516e-06 



-7.7554691781e-02 

3.76l8465057e-04 

-4.334 5626096e-07 

q 4 


2.4783745998e-03 

-1.21 23 l 65585e-05 

1 .-1 181424491 e— 08 




A' 7/4 H SO 4 ( n=5) 


Qo 


9.9601477451e-01 

-1. 1296951 29 le-02 

1.7452716744e-05 



1.4845831451e+00 

-9.3756384679e-03 

1.8989492030e-05 

a 2 


-2.9273249432e-01 

1.3469653174e-03 

-2.27189 17359e-06 

a 3 


2.1988585949e-02' 

-8.14 77260064 e-05 

I.0470356871e-07 

a A 


-7.53971 12432e-04 

2.2096863963e-06 

-1.6051821078e-09 

a 5 


9.7902066682e-06 

-2. 1560396221e-08 

-3.3248336670e-12 

a In7 = 

EIU 

ait(T)m k/2 , Read 1.1550607531e+01 as 1.1550607531x10'. All 


rr cn / u r\ — * — 4 - ww “ lu “ dttivuy tor cne 

rt>bU.,/h_>U system, was not parameterized but instead tabulated in the computer 
codes. r 


.\H 3(g) + H 


(aq) 



(13) 


where Kj is the equilibrium constant and (g) and (aq) 
refer to gas- and aqueous-phase species, respectively. 
The first dissociation step of H_»SO.| tH_»SO., t |) — > 
HSO. t(aq , + is assumed to be complete. 

Temperature-dependent water and solute activity co- 
efficients in EQUISOLV II are modified at cold tem- 
peratures (190°K-260°K) by parameterizing data from 
the model of Clegg et al. [1998a] The temperature- 
dependent water activities of five electrolytes in\olved 
in equations (11)-(13), (H ,SO,. HXO,.' (XH,),SO.,. 
XH.iXOj. XH , HSO.|) are improved in EQUISOLV II. 
For consistency with previous work by Tabazadeh et al. 
[1997a, 1997b] the binary solution compositions are pa- 
rameterized into 


m = y,(a w ) 


(T - 190) [yj(a w ) - y t utw )1 
70 


(14) 


" here m is the molality of the binary electrolvtes. a w is 
the water activity (relative humidity expressed as frac- 
tion), and y is a function of a vv only (see Table 5 ). 
Mean binary activity coefficients ( 7 ? ) were parame- 
terized by converting mole fraction activity coefficients 
(0 of Clegg et al .' s [1998a] model into molality base 
( 7 ) values according to equation (13), and using = 

(" 1 , + 7 [Pitzer, 1991]: 



where i/ + and 1 /_ are the stoichiometric coefficients of 
the binary electrolytes (for example. u + = 2 and = 1 
for (XH,),,SO,). M w is the molecular weight of water, 
and summation is over all solute species. The mean' 
activity coefficients were then fitted into simple poly- 
nomial functions (given in Table 6 ) of temperature and 
molality (where a*, and 3 k are polynomial coefficients): 
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Cases 

Background 

Volcanic 

Background 

Polluted 


Background 


Conditions 11 

nh 3 

NH 3 , ppt 

H 2 SO 4 , ppt 

hnc> 3 , ppt 

H 2 O, ppm 

Illustration 

H 2 S 04 









Stratospheric Simulation {Ptotal— $0 mbar ) 




r = 0.0 

0.0 

500.0 

10 , 000.0 


5.0 


Figure la and lb 

r = 0.0 

0.0 

20 , 000.0 

10 , 000.0 


5.0 


Figure 2a and 2b 


Upper Tropospheric Simulation (Ptotal -200 mbar ) 




r = 0.0 

0.0 

100.0 

100.0 

5.0, 

50.0, 

500.0 

Figure 3a 

r = 1.0 

100.0 

100.0 

100.0 

5.0, 

50.0, 

500.0 

Figure 3b 

0 

c-i 

11 

200.0 

100.0 

100.0 

5.0, 

50.0, 

500.0 

Figure 3c 

r = 0.0 

0.0 

100.0 

2 , 000.0 

5.0, 

50.0, 

500.0 

Figure 4a 

r = 1.0 

100.0 

100.0 

2 , 000.0 

5.0, 

50.0, 

500.0 

Figure 4b 

r = 2.0 

200.0 

100.0 

2 , 000.0 

5.0, 

50.0, 

500.0 

Figure 4c 


HNOz Vp take (Ptotal =200 mbar ) 





r = 0.0 - 2.0 

0.0 - 200.0 

100.0 

100.0 

5.0, 

50.0, 

500.0 

Figure 5 

r = OC 

100.0 

0.0 

200.0 

5.0, 

50.0, 

500.0 



a Ptotal is total atmospheric 
Laaksonen et a/., 1997). 


pressure, and the mixing ratios used are atmospheric observations [Tabazadeh et at., 1998; 


ln 7 -° = X>k( T )m k/2 , 


k = 0 


o k (T) = A + AT+AT 2 . 


The above parameterizations induce relative errors 
U6a) of no more than a few percents for both water and so- 
lute activities. The largest errors occur at low molal- 
(16b) ity regions where ion activities usually exhibit a com- 



(b) aerosol volume 


7 }Z 7 5 1 98 37 28 28 94 ’9 52 1 4 39 '0 62 7 90 5.91 

Correspond riQ Rh { % ) 


O O O w '* AiM2 ( Ce 99 el 0l * 1998 °) 

□ □ □ w2, At M2 (Ciegq et ol., 1998a) 

O O O vol “ AIW2 C c ’ e< 39 et ° 1 ' 199S °) 

■ this *ork (aPCm) 



_ — . — - this work (mod' f >ec CQutSOtv d 
ToDazcdeb es ol (199*) 
Carstaw el o' ( 1 995b) 

Luo et oh (1995) 


^ — o : 


186 188 190 192 194 ’96 '98 200 202 

Temperature (K) 


73.07 5198 37 28 26 94 19 62 1 4 39 ’0 62 7 9C 5 91 

Corresponding Rrt (5»; 


Figure 1. Comparison of various model predictions of (a) weight percents and (b) aerosol 
volume under background conditions in the stratosphere (see Table 7 for model conditions). 
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(b) aerosol volume 


6198 37.28 26.34 19 52 *4 39 1 0-62 7 90 5.9! 

Co'respond>3 Rh {%) 



73.07 51 98 37 28 26 9* ’9 62 14.39 t 0 62 7 90 59' 

Corresponding Rh { 7 .) 


Figure 2. Comparison of various model predictions of (a) weight percents and (b) aerosol 
volume under volcanic conditions in the stratosphere (see Table 7 for model conditions). 


plicated behavior. In addition to the improvement of 
binary water and solute activity coefficients, equilib- 
rium constants of equations (1 1)-( 13) were obtained 
from Clegg et al [1998a] and converted into suit- 
able molality-based units for use in EQUISOLV II. 
These modifications ensure perfect agreement of EQUI- 
SOLV II and AIM’2 at a binary level where electrolyte 
solutions are made ofHoO plus only one electrolyte. 

4. Model Intercomparison and 
Evaluation 

4 . 1 . h 2 so 4 /hno 3 /h 2 o 

For the stratospheric system of H^SO^HNOa/I-UO, 
the compositions obtained from APCM are compared 
with previous formulations [Carslaw el al, 1995b; Luo 
et al 1995: Tabazadeh et al.. 1994] under background 
and volcanic states (see Table 7 for model conditions 
used). Results from model imercomparisons between 
six different ternary aerosol models are shown in Fig- 
ures 1 and 2. For the background stratosphere, APCM 
with Diff correction terms included is in good agree- 
ment with AIM2 [Carslaw et al. } 1995a]. The predicted 
weight percents by APCM overlap with AIM2 almost 
exactly throughout the entire temperature range shown 
in Figure la. whereas previous models deviate from 


AIM2 either at the initial, in the middle, or at the 
final stage of aerosol growth (Figure lb). Under vol- 
canic conditions (H2SO4 increases from 0.5 to 20ppb). 
all models produce nearly identical results with only 
slight variations (Figures 2a and 2b). 

4 . 2 . NH3/H2SO4/HNO3/H2O 

In the upper troposphere, aerosol compositions are 
examined for both background and polluted states (see 
Table 7 for model conditions). Under polluted con- 
ditions, HN0 3 mixing ratios are elevated well beyond 
background values of lOOppt [Laaksonen et al., 1997], 
mainly because of convective transport of polluted boun- 
dary layer air directly into the upper troposphere. Vari- 
ations in the aerosol composition predicted by three 
models, APCM, modified EQUISOLV II (with the im- 
proved activity data), and AIM2 [Clegg et al. , 1998a], 
are illustrated in Figures 3 and 4. The points of ice sat- 
uration are marked as vertical dotted lines in Figures 3 
and 4. For all cases studied and compared, APCM 
yields compositions that are in close agreement with 
A I M2. The modified EQUISOLV' II with new activ- 
ity data performs better in the regions where relative 
humidity (RH) is high. The deviations at lower rela- 
tive humidity are likely caused by the uncertainties in 
mixed activity coefficients calculated from simple mix- 
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ing rules [Bromley, 1973; Kusik and Meissner , 1978] in 

EQUISOLV II. 

We have also compared predicted molalities of indi- 
vidual components, NH*, H2SO4 (= HSO^ + SOlJ ). 
and NOJ . in solution (instead of weight percents of 
wj and wj) and found that the differences between 
APCM and AIM2 are no more than a few percent ex- 
cept at the very low or high solute concentration re- 


gions. In other words, predictions by APCM are less 
accurate only in regions where weight percents approach 
the lower or upper limits (1% and 85%). Nevertheless, 
for very low weight percents, concentrations are too di- 
lute to be significant, whereas for very high weight per- 
cents, relative humidities are sufficiently low that am* 
moniated salts, such as (NH 4 )_>S0 4 or (NH 4 )3H(S0 4 )„> 
(letovicite), would most likely precipitate in solution 
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Figure 4. Variation of aerosol compositions as a function of temperature (or relative humidity) 
under polluted conditions in the upper troposphere (see Table 7 for model conditions), (a) 
r =0.0. (b) r = 1.0, and (c) r = 2.0. 


[Tahazadeh and Toon , 1998]. For example, crystalliza- 
tion of (NH.,)- 2 S 04 occurs at - 35% RH at room tem- 
perature [Xu at «/.. 1998]. At colder temperatures, crys- 
tallization of (NH 4 )->SO, will probably occur at slightly 
higher RHs. In addition, for very concentrated solu- 
tions (w t + w» > 85%), aerosol compositions predicted 
by AIM2 are often outside the range of model validation 


[Clegg et a/., 1998a]. For the reasons mentioned above, 
the limits imposed on the weight percents in the APCM 
are not a serious drawback. Also, as shown in Figures 3 
and 4, the APCM compositions are most accurate near 
regions where ice reaches saturation in the atmosphere. 

In addition to aerosol composition the APCM can be 
used to calculate the extent of HNO3 uptake by up- 
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Figure 57 Uptake of HN0 3 in sulfate-based aerosols as a function of sulfate neutralization 
under background conditions in the upper troposphere (see Table 7 for model conditions), (a) 
H.,0 = oppm, (b) H a O = 50ppm, and (c) HjO = SOOppm. (Dashed lines represent the uptake 
by pure NH4NO3 aerosols (i.e.. H2SO4 = 0)) 


per tropospheric aerosols. Aircraft field experiments 
have shown that ammoniated particles are abundant 
in the upper troposphere [Talbot et al, 1996, 199S: 
Tabazadeh et al., 1998). Here we examine how HNO3 
uptake may be affected by the presence of ammoni- 
ated aerosols in the upper troposphere. For simulations 
the HjSO, and HNO3 mixing ratios are fixed at back- 
ground levels of 10C ppt each. NH 3 concentrations in 
the APCM are varied by increasing (or decreasing) r 
to account for changes in solution neutrality (or pH), 
ranging from pure HjSO.i solution droplets (r = 0.0) to 
fully ammoniated systems consisting of only (NH^oSO^ 
(r = 2.0) isee Table 7 for model conditions used). 


In Figure 5, the uptake of HNO3, expressed as a frac- 
tion (i.e., ratio of concentration in the liquid phase to 
the total initial concentration), is shown for three differ- 
ent assumed water vapor pressure profiles. As expected, 
partitioning of HNO3 in sulfate-based aerosols depends 
strongly on solution neutrality. In general, greater up- 
take occurs at lower temperatures and higher relative 
humidities. For example, at temperatures lower than 
210°K (which corresponds to a relative humidity of 78% 
for P tota i = 200mb and H 2 0 = 50ppm), a significant 
fraction (> 90%) of HNO3 resides in fully neutralized 
ammoniated solutions, compared to only < 20% in the 
pure sulfate system (Figure 5b). It has been shown that 
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high levels of HNO 3 in solution may cause precipita- 
tion of arnmoniated and/or nitrated salts, which could 
change the mode of ice formation from homogeneous to 
heterogeneous nucleation [ Tabazadeh and Toon , 1998]. 

5. Conclusions 

The ion interaction model of Clegg et ai [1998a] for 
the system of NHa/HoSO^/HNOa/H^O has been pa- 
rameterized into a compact model well suited for in- 
corporation into large-scale atmospheric models. The 
aerosol physical chemistry model (APCM) reproduces 
the AIM‘2 results of Clegg et al. [1998a] for a wide 
range of conditions in the upper troposphere and lower 
stratosphere. Model intercomparisons show that for the 
ternary system of H2SO4/HNO3/H0O solution compo- 
sitions obtained from APCM are in better agreement 
with those obtained from AIM2 than previous formula- 
tions. For the quaternary system of NH3/H2SO4/HNO3 
/H 2 O, APCM results are also in good agreement with 
AIM2 predictions, particularly near the regions of ice 
saturation where the influence of arnmoniated particles 
on the ice nucleation process is of interest. Extension 
of APCM to include other features such as calculations 
of deliquescence relative humidity and precipitation of 
solids in solution are under way. 
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Abstract. 

Satellite observations of water vapor and aerosol extinction along with temperature 
trajectory calculations are analyzed for the Southern Hemisphere winter of 1992 
in order to determine the onset, extent and duration of dehydration within the 
polar vortex. Our investigation utilizes measurements of water vapor from the 
Microwave Limb Sounder (MLS) and extinction from the Cryogenic Limb Array 
Etalon Spectrometer (CLAES) both onboard the Upper Atmosphere Research 
Satellite (UARS). Evidence for dehydration is seen on potential temperature 
surfaces from 420-520 K (approximately 16 to 22 km). The horizontal extent of 
the dehydrated area at 465 K encompasses up to 35% of the total vortex area. 
Based on MLS observations and temperature statistics, the onset of dehydration 
occurs between late June and early July, and dehydrated regions persist well into 
November. A comparison of CLAES aerosol extinction measurements and model 
calculations of extinction suggests an average ice particle number concentration 
and size of 10~ 1 2 -10 -3 cm -3 and 10-30 ^m, respectively. We show that the 
difference in timing of the onset of dehydration found here an in a recent analysis 
of Polar Ozone and Aerosol Measurement III (POAM) observations is due to the 
latitudinal sampling pattern of the POAM instrument. 


1. Introduction 

Polar stratospheric clouds (PSCs) form in the low temper- 
atures of the winter polar vortex and have a direct influence 
on atmospheric chemistry by their involvement in heteroge- 
neous ozone chemistry {Solomon, 1999]. PSCs also act as 
a sink for gas phase water vapor and nitric acid resulting in 
stratospheric dehydration [e.g., Kelly el al . , 1989; Hofmann 


and Deshler , 1991; Nedoluha et al . , 2000] and denitrifica- 
tion [e.g., Toon et al., 1986; Santee et al, 1999; Tabazadeh 
et al., 2000a]. 

Understanding the physical properties of the clouds and 
the meteorological conditions under which they form is im- 
portant for studies of the polar vortex. Particle size and com- 
position influence the optical properties of the clouds, which 
affect remote sensing [Toon et al., 1990; Browell et al. 
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1990]; composition of the aerosols impacts constituent bud* 
gets and chemical reactions [ Solomon , 1999]; and chemical 
processes which are dependent upon relative humidity will 
be affected by the deficit of water vapor due to dehydration 
[Hanson et al , 1994]. PSCs are separated into two types 
[WMO, 1999]. Type II PSCs are composed of water ice and 
form at temperatures below the ice frost point. Type 1 clouds 
form at temperatures above the ice frost point and consist of 
both liquid and solid nitric acid*containing cloud particles. 
We focus our study on the atmospheric dehydration process 
which occurs as a result of ice particle growth and sedimen- 
tation, thereby removing water from the atmosphere. 

The Airborne Antarctic Ozone Experiment (AAOE) in 
August-September 1987 show ed the first evidence of strato- 
spheric dehydration in the Antarctic vortex [Kelly et al, 
1989], Balloon measurements in the Antarctic have also 
shown profiles with water vapor depletion [Hofmann et al, 
1991 \ Rosen et al , 1988; Votnel et al, 1995]. Using balloon 
measurements over McMurdo in 1994, Voniel et al [1995] 
conclude that dehydration starts around mid June and per- 
sists into November. Overall, only a limited number of in 
situ measurements have been made over the Antarctic re- 
gion over the course of several winters. Satellite observa- 
tions can provide a more extensive view of the polar re- 
gions. Springtime dehydration has been noted in the analysis 
of Halogen Occul-tation Experiment (HALOE) data [Pierce 
et al, 1994; Rosenlof et al , 1997]. Santee et al. [1995] re- 
port on the interhemispheric differences in wintertime water 
vapor from the Microwave Limb Sounder (MLS) and Ricaud 
et al [1995] deduce PSCs from MLS water vapor and Cryo- 
genic Limb Array Etalon Spectrometer (CLAES) extinc- 
tion measurements from several days in August-September 
1992. Evidence for wintertime dehydration is also found in 
measurements from the Atmospheric Trace Molecule Spec- 
troscopy (ATMOS) instrument [Manney et al, 1999], and 
Nedoluha et al [2000] report on Polar Ozone and Aerosol 
Measurement III (POAM) observations which show dehy- 
dration occurring in the Southern Hemisphere (SH) winter 
of 1998. 

In our study, we examine satellite measurementsof strato- 
spheric water vapor and extinction in the SH winter polar 
regions to identify the onset, extent and duration of dehydra- 
tion. We correlate water vapor concentrations with aerosol 
extinction and use temperature trajectories to examine where 
and when the dehydration process occurs. We compare mea- 
sured and calculated extinction to infer ice particle size and 
number density in the dehydration column. 


2. Observations 

Water vapor data for this study come from the MLS in- 
strument onboard the Upper Atmosphere Research Satellite 
(UARS). We use the prototype water vapor product (version 
104) described by Pumphrey [1999], which gives water va- 
por mixing ratios in the stratosphere from 0.1 to 100 hPa 
with a vertical resolution of 3-4 km. Although this prod- 
uct is not an official UARS project MLS product, its validity 
and usefulness in scientific studies has been demonstrated 
[Pumphrey et al, 2000; Pumphrey, 1999; Manney et al, 
1998]. The water vapor measurements are available from 
September 1991 to April 1993, which includes only one SH 
winter. 

The CLAES instrument, also onboard UARS, measures 
aerosol extinction coefficients [Mergenthaler et al, 1997] 
which have been used to identify PSCs [Mergeruhaler et al, 
1997; Ricaud et al, 1995; Massie et al, 1994] and to de- 
termine volume and surface area densities of type I PSCs 
[Massie et al, 1997, 1998]. In identifying PSCs, extinction 
values are generally used qualitatively and threshold values 
for PSCs, depending on pressure, are determined relative to 
temperature and pre-winter observations. No information on 
particle size and concentration is given with the aerosol ex- 
tinction data product. CLAES was operational from Septem- 
ber 1991 until May 1993. Our analysis uses the CLAES ver- 
sion 8 extinction coefficients at 780 cm” 1 . 

MLS and CLAES measurements are nearly coincident 
in time and space, each with a horizontal resolution of 
^400 km. The vertical profiles of water vapor and extinc- 
tion have been interpolated from the standard UARS grid to 
potential temperature surfaces using United Kingdom Mete- 
orological Office (UKMO) temperatures. UKMO temper- 
atures are also used for the analysis in sections 2 and 3. 
Comparisons of the UKMO temperatures with radiosonde 
observations show that systematic biases are less than 1 K 
throughout the SH winter [Manney et al, 1996]. 

We focus on the UARS SH viewing periods for one year, 
starting in April 1992 ( 1-30 April 1992, 2 June- 12 July 1992, 
14 August-20 September 1992, 30 October-28 November 
1992, 10 January-8 February 1993, and 20 March- 1 April 
1993). We note that there are significant data gaps in the 
June-July viewing period due to problems with the UARS 
solar array which required the MLS and CLAES instruments 
to be turned off. Data availability limits this study to only 
one SH winter and we therefore do not assess year-to-year 
variability of dehydration. Although the SH wdnter of 1992 
was not necessarily a typical year since there were high lev- 
els of stratospheric aerosol from Pinatubo, the duration of 
cold temperatures and minimum temperatures were typical 
of wintertime conditions. 
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The water vapor annual cycles on potential temperature 
surfaces, compiled from the high-latitude MLS data, are 
shown in Figure 1. Daily averages are from data in the 
70°S-80°S latitude range for one year starting in April 1992. 
The 420 K potential temperature level is below 100 hPa in 
January through March, resulting in a lack of observations 
on that surface during that time. Potential temperature sur- 
faces from 420 to 585 K show a minimum in water vapor 
in the winter months. These potential temperature surfaces 
cover pressure levels of 100 to 15 hPa. Since the water vapor 
data set does not extend below 100 hPa, we cannot observe 
the lower altitude limit of the dehydration. Nedoluha et ai 
[2000] show evidence for dehydration down to 12 km. At 
the higher potential temperature surfaces, the maximum in 
the annual cycle is found in late winter (August-September), 
The depleted winter values at and below 585 K are a result 
of cold temperatures found within the polar vortex. Average 
temperatures in July on the 465 K potential temperature sur- 
face are near 185 K, below the ice frost point (188 K). Ice 
clouds form at these cold temperatures and deplete the gas 
phase water. 

The daily averages of aerosol extinction from CLAES are 
also shown in Figure 1. In general, they are anti-correlated 
with the water vapor in the winter months from 420 to 585 K. 
High extinction values are seen from mid June through late 
August when the water vapor mixing ratios are low, rein- 
forcing the observation that the deficit in gas phase water 
coincides with the presence of PSCs. Extinction values in 
September are lower than the pre-winter values. This fea- 
ture in the annual development of extinction has been seen 
by other observations and is due to the cleansing of the at- 
mosphere by the subsidence of PSCs and by descent in the 
vortex [Kent et al.\ 1985; Thomason and Poole , 1993; Ran- 
dall et ai , 1996; Nedoluha et ai , 2000] 

Due to limited MLS water vapor measurements in late 
June and early July, we eannot observe the precise onset of 
dehydration. The measurements indicate that dehydration 
begins sometime after June 1 8. Extinction values are already 
high prior to June 1 8 due to Type I PSCs [Tabazadeh et ai , 
2000a). which form at temperatures higher than the ice frost 
point. A deficit in the water vapor mixing ratio is clearly 
seen in the July 10-12 measurements for the 420 to 585 K 
surfaces. At 465 K, daily average values have dropped from 
4 ppmv in mid June to 2.7 ppmv on July 10, indicating that 
the onset of substantial dehydration occurred sometime af- 
ter mid June and before 10 July. Water vapor values at the 
420 K and 465 K levels drop further between mid July and 
early August when the next observation period begins. In 
August-September, daily mean values are between 2.3 and 
2.8 ppmv at 465 K. Extinction values are still elevated in mid 
August but have diminished by early September. Measure- 


ments on the 585 K surface show a drop in the water vapor 
and an elevation of extinction in mid July, but values have 
recovered by August, indicating that the initial drop in water 
vapor seen at this level was reversible or that mixing ratios 
are replenished by the descent of moister air from above. 

Figure 1 depicts only daily averages in the 70°S-80°S lat- 
itude range. Additional information can be found in the dis- 
tributions of water vapor and extinction measurements. Fig- 
ures 2a and 2b show distributions of water vapor and aerosol 
extinction in the region of 50°S-80°S for different temper- 
ature bins (T < 1 85 K, 1 85 <T <188 K, 1 88 <T< 195 K, and 
1 95 <T <210 K (points with temperatures greater than 2 10 K 
are excluded). We also show the distributions of water vapor 
and extinction for the 70°S-80°S latitude bin (for all temper- 
atures). TheT< 185 K bin represents air below the ice nucle- 
ation point (typically 3-4 K below the frost point [Tabazadeh 
et ai , 2000b; Chang et ai , 1999]) and the 1 85 <T< 188 K 
bin represents observations between the frost point and the 
nucleation point. The distributions for data taken between 
April and November 1992 are separated according to the 
UARS yaw cycle. 

We use the high latitude April measurements, represent- 
ing pre-wintertime observations, to define low water vapor 
and high extinction values later in the paper. The water vapor 
distribution peaks at 4.3 ppmv in April and all the extinction 
values in this month are less than 0.0009 km _; . In June- 
July, the mode of the water vapor distribution for T< 185 K 
is 2.3 ppmv and the i 85 <T< 188 K observations have two 
peaks at 2.3 ppmv and near 4 ppmv. For T> 1 88 K the distri- 
butions are centered at 4 ppmv and fewer than 10% of the ob- 
servations are less than 3.5 ppmv, indicating that by this time 
the dehydrated area has not yet spread beyond the stability 
region of ice particles. Therefore, we define dehydrated ar- 
eas to have water vapor mixing ratios less than 3.0 ppmv 
based on the April measurements, and “high” aerosol ex- 
tinction is defined to be values larger than 0.001 krrT l . The 
aerosol extinction threshold is slightly higher than the value 
used to identify PSCs at 46 hPa (0.00075 km“ ! ) by Mergen - 
thaler et ai [1997]. Extinction values over 0.002 km“ l have 
large uncertainties due to the high optical depth of the PSCs 
along the limb viewing path [Mergenthaler et ai , 1997], 
However, the uncertainties do not affect the classification of 
“'high” aerosol amount. The June-July distributions show 
that nearly all the observations which coincide with temper- 
atures below 185 K are dehydrated. 25% of the observa- 
tions with I85<T< 188 K are dehydrated. For temperatures 
above 188 K, there is no evidence of significant depletion of 
the water vapor. 

The extinction distributions for the June-July observation 
period show evidence of both Type I and Type II PSCs. 
The temperature bins with T < 1 95 K all show observa- 
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tions of extinction above the April distributions. For the 
1 85 <T< 188 K temperature bin, 75% of the observations 
have extinctions larger than the pre-winter values. This is 
significantly larger than the 25% of water vapor measure- 
ments that show a deficit in this temperature bin, suggesting 
that type I PSCs are contributing to the elevated extinction. 
The coldest temperatures coincide with the highest extinc- 
tion values. All of the observations with T> 195 K show low 
extinction values, less than 0.001 km" 1 . 

By the August-September viewing period, the distribu- 
tions of water vapor mixing ratios for all temperature bins 
with T<195 K peak near 2.2 ppmv. For 195<T<2iOK, the 
majority of the observations are centered around 4 ppmv but 
there are measurements with mixing ratios less than 3 ppmv 
in this temperature range as well. The extinction distribu- 
tions during this time have a small percentage of observa- 
tions with extinctions above 0.001 km* -1 and these coincide 
with the coldest temperatures (T< 188 K). Thus by this time 
in the winter, there exist dehydrated mixing ratios that do 
not directly coincide with low temperatures or high extinc- 
tion values. 

Although temperatures at high latitudes are warmer by 
November, water vapor mixing ratios have recovered only 
slightly from the dehydration. The distribution for 70°S- 
80°S peaks at 3.0 ppmv. While some water vapor recovery 
has taken place in November, the overall mixing ratios are 
still substantially below the April mode value (~4.3 ppmv) 
before any dehydration occurred. 

3, Extent and Timing of Dehydration 

In this section we quantify the horizontal extent of the 
dehydration, its development throughout the winter months, 
and its relationship to the areas of cold temperatures, high 
extinction, and the vortex. To do this, we calculate the area 
encompassing the vortex and the regions of dehydrated air, 
high extinction, and cold temperatures. Criteria for dehydra- 
tion (mixing ratios less than 3.0 ppmv) and “high” extinction 
(larger than 0.001 km"* 1 ) are determined from the frequency 
distributions in the early winter months and before the onset 
of winter (Figures 2a and 2b). The water vapor distributions 
for early winter show two peaks, separated near 3 ppmv, and 
April measurements all exceed 3 ppmv (Figure 2a). The ex- 
tinction distributions in early winter show- a distinct separa- 
tion at 0.001 km" 1 and pre-winter values are all below that 
threshold (Figure 2b). The cold temperature criteria are the 
frost point of ice (T < 1 88 K) and the ice nucleation tem- 
perature (typically 3 K less than the frost point. T< 1 85 K 
[Tabazadeh et ai , 2000b; Chang et at. , 1999]). We use a 
scaled potential vorticity value of 1.4 s ~ 1 to approximate the 
vortex boundary. Figure 3 shows the calculated daily areas 


from April through November on the 465 K potential tem- 
perature surface. 

The area of temperatures below the frost point of ice be- 
gins to develop in mid June. At this time high extinction 
values are present but there is no evidence of dehydration . 
The area of high extinction values drops on June 24 and then 
increases again. Tabazadeh et al. [2000a] show that deni- 
trification at 450 K occurred over 12 days from June 12-24. 
Temperatures below the ice nucleation point are present over 
a small area in mid June but do not persist or develop over 
larger areas until after June 27. By July 10, the area of dehy- 
dration coincides roughly with the area of temperatures be- 
low the ice nucleation point and is similar to the area of high 
extinction. When water vapor measurements are available 
again in mid August, the dehydrated area exceeds the cold 
temperature area and persists at 12 x 10 6 km 2 through most 
of September, consistent with observations by [Vdmel et ai % 
1995] in which they find that the water vapor profiles from 
early August and early October in 1994 at McMurdo Station 
are similar. From mid August through late September the de- 
hydrated area is roughly 30%-35% of the total vortex area. 
Cold temperature areas decrease rapidly at the end of Au- 
gust. The area of high extinction is much less than the dehy- 
drated area during August and September, indicating that the 
dehydration persists after the ice clouds have precipitated. 
The dehydrated areas persist and remain at nearly a constant 
value (from early August to late September) after the area of 
high extinction (indicating the presence of clouds) has van- 
ished in early September and after the vortex area begins to 
recede in early September. The dehydrated area decreases 
rapidly after November 15, and the water vapor values start 
a rapid recovery. The water vapor mixing ratios at 465 K 
(Figure 1) show an increase between mid September and the 
beginning of November similar to what is observed for nitric 
acid [Santee et ai , 1999]. The dehydrated area in Figure 3 
does not show this recovery occurring since the area will not 
start to decrease until mixing ratios surpass 3 ppmv, indicat- 
ing that mixing ratios are attaining values similar to those 
present before temperatures fell below the ice frost point. 

4. Temperature Statistics 

The interpretation of the time development of the area 
of dehydration compared to the area of cold temperatures 
(Figure 3) is limited in that it does not include the temper- 
ature history of air parcels. Tabazadeh et al. [2000a] uti- 
lize the concept of “PSC lifetimes” to examine how long an 
air mass persists at temperatures less than the condensation 
points of nitric acid trihydrate (NAT), nitric acid dihydrate 
(NAD), and water ice. In this section, we use temperature 
statistics from National Centers for Environmental Predic- 
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lion (NCEP) trajectory cumulations to look at the onset of 
temperatures which meet the ice nucleation conditions for 
the formation and growth of water ice clouds. 

For each day in June through September, 20 points were 
homogeneously distributed within the 188 K temperature 
contour at 50 hPa. Twenty forward and backward trajec- 
tories were run isentropically from the initialized points for 
14 days. Diabatic descent, which can impact the Antarc- 
tic lower stratosphere over a 14-day period in midwinter, is 
neglected in these calculations. The statistics of the trajec- 
tory analysis are shown in Figure 4. The solid line depicts 
the percentage of trajectories for each day that have passed 
through the ice nucleation criterion (reaching temperatures 
of about 185 K or below). Symbols show' the average, min- 
imum and maximum time that the ensemble of trajectories 
spends below the ice frost point (188 K) after the nucleation 
criterion has been met. The temperature statistics suggest 
that the onset of dehydration occurred at the end of June 
when a significant number of the trajectories met the nucle- 
ation criterion for the formation of ice particles. By early 
July, 60 to 90% of the parcels encountered temperatures be- 
low 185 K, suggesting that by this time ice particles nucle- 
ated in nearly all the air parcels initialized within the 188 K 
temperature contour. In addition, the lifetime of ice parti- 
cles nucleated was long enough (5 days or longer) to provide 
sufficient time for growth and sedimentation of particles to 
lower altitudes. 

5. Model Extinctions 

Next we focus on the extinction measurements from 
CLAES. Previous studies have used these measurements 
to identify PSCs [Mergenthaler et al, 1997; Ricaud et al, 
1995; Massie et al, 1994]. Here we use Mie calculations to 
compute aerosol extinction due to ice particles. The calcu- 
lated extinction values are compared to the CLAES observa- 
tions to infer ice particle size and number density. 

Our calculated extinction depend upon temperatures, pres- 
sure, water vapor mixing ratio and fraction of aerosols nucle- 
ated into ice particles. To compute the extinction, we have 
distributed the volume of condensed water (as ice) at equi- 
librium over a fixed and assumed number density of ice par- 
ticles. The size of the ice particles obtained in this manner, 
for a given number density, were used as input parameters 
into a Mie code to convert the volume and size information 
into extinction units [Toon et al, 1990; Massie et al , 1994]. 
The CLAES observations are compared with extinction cal- 
culated using different assumptions as to the number density 
of ice particles based on the range of number densities ob- 
served in the in situ data [Hofmann et al , 1991 ]. Maps of the 
calculated extinction are formed by computing the extinc- 


tion for the temperatures at the geographic locations of the 
CLAES measurements when the temperature is below the 
threshold for type II PSCs (188 K). At locations where this 
temperature criterion is not met, the observed CLAES val- 
ues are used in constructing the extinction maps. The maps 
depict what the CLAES instrument would have observed if 
type II PSCs were present in the locations where they are 
predicted to form and if the instrument measured perfectly. 
We examine days from early in the winter since once con- 
densation and ice particle formation begins, the amount of 
water vapor available changes. For the early wintertime 
conditions, we use a mixing ratio of 4.1 ppmv, based on 
the mean values in Figure 1. Figure 5 shows maps of the 
CLAES data and calculated extinction for three different as- 
sumptions of ice particle density. The figure shows that an 
ice particle number density in the range of 10“ 2 -10~ 3 cm“ 3 
best matches the CLAES observations. The comparison be- 
tween the observed and calculated extinction is good despite 
the uncertainty in the CLAES measurements for high val- 
ues of extinction. Based on the modeling analyses of Jensen 
andToon [1994], the low ice number density implies that the 
clouds were formed in regions with very slow cooling. 

The size of ice particles associated with the concentra- 
tions of 10"MCr 3 crrT 3 is shown in Figure 6. The distri- 
butions are formed from the calculated sizes for the days in 
Figure 5. Ice particle sizes are on the order of 10-30 fim. 
These should be considered general values since the real at- 
mosphere is likely to contain a distribution of ice particle 
number densities and sizes. In general our inferred ice par- 
ticle number density and size are in good agreement with 
balloon observations over McMurdo station [Hofmann et al , 
1991]. 

6. Comparison of dehydration from POAM 
and MLS 

A recent study of water vapor and extinction data from 
POAM III examines dehydration in the Antarctic for the SH 
winter of 1998 [Nedoluha et al 2000]. POAM is a solar oc- 
cultation instrument and makes 14 measurements per day in 
each hemisphere around a latitude circle. The latitude sam- 
pled varies semi-annually as shown in Figure 7. The POAM 
instrument measures water vapor and extinction with a high 
vertical resolution (~1 km) in the lower stratosphere, mak- 
ing it useful for studies of dehydration and PSC occurrence. 
Nedoluha et al [2000] report dehydration occurring over a 
6 week period starting in mid July with the drop in water 
vapor coincident with the timing of the minimum tempera- 
tures dropping below the ice frost point in air parcels sam- 
pled by the POAM instrument and with the occurrence of 
a significant number of extinction measurements indicating 
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the presence of type II PSCs. 

While the data used in our study cover a different calen- 
dar year, we note that our analysis of the water vapor mix- 
ing ratios from MLS and temperature statistics suggests that 
the onset of dehydration begins prior to early July, and it is 
therefore necessary that we relate our conclusions to those 
of Nedoluha et al [2000]. MLS mixing ratios from 70°S to 
80°S drop significantly (from 4 to 2.7 ppmv) from mid June 
to July 10 (Figure 1 ). The area of dehydration and cold tem- 
peratures begins to develop in late June (^Figure 3). Vomel 
ei al. [1995] suggest that dehydration and the cold temper- 
ature regions develop initially in the interior of the vortex 
and expand outward. Thus the first observation of dehy- 
drated air depends upon the latitude being sampled. From 
late June through mid July, the POAM instrument observes 
latitudes from 65°S-70°S. The MLS water vapor measure- 
ments in this latitude range do not show evidence of signifi- 
cant dehydration in 10-12 July (Figure 7), whereas substan- 
tial dehydration is evident in the MLS data closer to the pole 
(70°S-80°S) for those dates. Because MLS shows no signs 
of dehydration in mid June and substantial dehydration on 
10 July, we conclude that the onset of severe dehydration 
occurred sometime between mid June and 10 July. Further, 
temperature statistics provide additional support for the on- 
set of dehydration to be around late June. Therefore, if the 
meteorological conditions during the 1992 and 1998Antarc-‘ 
tic southern winters are comparable, then the onset of dehy- 
dration observed by POAM type instruments will always be 
offset as compared to the onset seen by MLS type instru- 
ments. By August, POAM is observing latitudes poleward 
of 70°S and deficits in the water vapor are apparent in the 
data during this period [Nedoluha ei al. . 2000]. Thus, while 
the analysis of the POAM data provides a much finer ver- 
tical picture and more continuous daily measurements than 
does the MLS data, it is important to note that the latitudinal 
sampling influences the determination of the onset of dehy- 
dration. 

7. Summary 

Examining satellite measurements of water vapor from 
the SH winter of 1992, we find evidence for dehydration due 
to the formation of type II PSCs on potential temperature 
surfaces from 420 to 520 K (16 to 22 km;. However, POAM 
observations of water vapor show that dehydration extends 
down to ~360 K (12 km). Thus the vertical extent of de- 
hydration is at least 4 km deeper than what MLS can see, 
according to the POAM data. The deficits in water vapor are 
seen in measurements from mid July persisting into Novem- 
ber. Mean mixing ratios decrease by approximately I to 2 
ppmv. Temperature statistics suggest that the onset of dehy- 


dration occurs in late June when a large number of cold air 
mass trajectories meet the ice nucleation criteria. Because 
the onset of dehydration initially occurs close to the pole, the 
latitudinal sampling of the atmosphere by different instru- 
ments must be taken into account for the determination of 
the dehydration onset. The development of dehydration over 
the course of the winter shows that initially (June-July obser- 
vations) dehydration is found only in areas which have cold 
temperatures (below 1 88 K). These areas coincide with those 
with high aerosol extinction. Area calculations show that 
regions which meet the cold temperature criteria and those 
which meet the criteria for dehydrated air expand through- 
out July. By August, the dehydrated regions are larger than 
those of the cold temperatures and high extinction. The de- 
hydrated mixing ratios persist after the elevated extinction 
and cold temperatures diminish. Results from a Mie code 
calculation in conjunction with the extinction measurements 
give a general estimate of aerosol properties. A number den- 
sity between I0" 2 - 10“ 3 cm" 3 best matches the observed ex- 
tinction and, our estimate of ice particle size is on the order 
of 10-30 pm. 
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